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. , Abstract. The chaotic or ordered character of orbits in galactic models is an important issue, since 

it can influence dynamical evolution. This distinction can be achieved with the help of the Smaller 
ALingment /ndex - {SALI). We describe here briefly this method and its advantages. Then we 
, apply it to that case of 2D and 3D barred galaxy potentials. In particular, we find the fraction of 

■ chaotic and ordered orbits in such potentials and present how this fraction changes when the 

7— I \ main parameters of the model are varied. For this, we consider models with different bar mass, bar 

thickness or pattern speed. Varying only one parameter at a time, we find that bars that are more 
massive, or thinner, or faster, have a larger fraction of chaotic orbits. 
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^ INTRODUCTION 

6 ■ 

The chaotic and ordered motion in dynamical systems is a very important issue in non 
- linear sciences in general and in dynamical astronomy in particular. The qualitative 
distinction of the motion is, generally, not that trivial and it becomes harder in systems 

. ^ | of many degrees of freedom. Thus, it is necessary to use fast and precise methods to 

^ ■ indicate the nature of the orbits. 

Many methods have been developed over the years trying to give an answer to this 
problem. The inspection of the successive intersections of an orbit with a Poincare sur- 
face of section (PSS) JH] has been used mainly for 2 dimensional (2D) maps and 2 
degrees of freedom Hamiltonian systems. One of the most common methods of chaos 
detection is the computation of the "Maximal Lyapunov Characteristic Number" (LCN) 
which can be applied for systems with many degrees of freedom. Another very 
efficient method is the "Frequency Map Analysis" [ 4] ,[ 5] . In recent years new methods 
have been introduced such as the study of spectra of "Short Time Lyapunov Characteris- 
tic Numbers" 1 6], |j] or "Stretching Numbers" [8], [9] and the "Spectral Distance" of such 
spectra ifioll . as well as the study of spectra of helicity and twist angles lfnll . lfl2ll . lfl3ll . In 



addition, Froeschle, Lega & Gonzi [ 141 introduced the "Fast Lyapunov Indicator" (FLI), 
while Vozikis, Varvoglis & Tsiganis [15] proposed a method based on the frequency 
analysis of "Stretching Numbers". 

Recently, a new, fast and easy to compute indicator of the chaotic or ordered nature 



of orbits, has been introduced Jl^l a nd was applied successfully in different dynamical 
systems 1 161. lfl7ll . lfl8ll . lfT9ll J20ll . ll2ill . ll22ll . ll23fl . the so-called "Smaller Alignment Index 
(SALI)", or, as elsewhere called Alignment Index (AI) [24], [25]. In the present paper, 
we first recall its definition and then we show its effectiveness in distinguishing between 
ordered and chaotic motion, by applying it to a barred potential of 2 and 3 degrees of 
freedom. 



DEFINITION OF THE SMALLER ALIGNMENT INDEX 

Let us consider the n - dimensional phase space of a conservative dynamical system, 
which could be a symplectic map or a Hamiltonian flow. We consider also an orbit in 
that space with initial condition P(0) = (jci(0),X2(0), ...,x n (0)) and a deviation vector 
v(0) = (dx\(0),dx2(0), ...,dx n (0)) from the initial point P(0). In order to compute the 
SALI for a given orbit one has to follow the time evolution of the orbit with initial 
condition itself and also of two deviation vectors vi (t) , vi{t) which initially point in two 
different directions. The evolution of these deviation vectors is given by the variational 
equations for a flow and by the tangent map for a discrete-time system. At every time step 
the two deviation vectors V\(t) and V2(t) are normalized and the SALI is then computed 
as: 



SALI(t) = min 



vi (0 + v 2 (t) 



vi (t) v 2 (0 



I vi (Oil I WO I 



(i) 



WOII IWO 

The properties of the time evolution of the SALI clearly distinguish between ordered 
and chaotic motion as follows: In the case of Hamiltonian flows or dimensional sym- 
plectic maps with n>2, the SALI fluctuates around a non-zero value for ordered orbits, 
while it tends to zero for chaotic orbits [ 16], [ 18], [ 19]. In the case of 2D maps the SALI 
tends to zero both for ordered and chaotic orbits, following however completely differ- 
ent time rates, which again allows us to distinguish between the two cases. In 2 and 3 
degrees of freedom Hamiltonian systems the distinction between ordered and chaotic 
motion is easy because the ordered motion occurs on a 2D or 4D torus respectively on 
which any initial deviation vector becomes almost tangent after a short transient period. 
In general, two different initial deviation vectors become tangent to different directions 
on the torus, producing different sequences of vectors, so that SALI does not tend to zero 
but fluctuates around positive values. On the other hand, for chaotic orbits, any two ini- 
tially different deviation vectors tend to coincide in the direction defined by the nearby 
unstable manifold and hence either coincides with each other, or become opposite. This 
means that the SALI tends to zero when the orbit is chaotic and to a non-zero value 
when the orbit is ordered. Thus, the completely different behavior of the SALI helps us 
to distinguish between ordered and chaotic motion in Hamiltonian systems with 2 and 3 
degrees of freedom and in general in dynamical systems of higher dimensionality. 



THE MODEL 



A 3D rotating bar model can be described by the Hamiltonian function: 

H=^(pl + Py+p 2 z )+V(x,y,z)-a b (xp y -ypx x ), (2) 

where x is the long axis, y the intermediate and z the short axis around which the bar 
rotates. The p x ,p y and p z are the canonically conjugate momenta. Finally, V is the 
potential, Qf, represents the pattern speed of the bar and H is the total energy of the 
system. 

The potential of the model V consists of three components: 
a) A disc, represented by a Miyamoto disc ll26ll : 

V„=- , GM ° (3) 



x 2 +y 2 + (A + Vz 2 +B 2 ) 2 

where Mo is the total mass of the disc, A and B are the horizontal and vertical scale- 
lengths, and G is the gravitational constant. 

b) The bulge is modeled by a Plummer sphere with the potential: 

y,=- /2 7\ , «> 

/x 2 +y 2 + z 2 + £ 2 



where e s is the scalelength of the bulge and Ms is its total mass, 
c) The triaxial Ferrers bar, the density p(x) of which is: 




where p c = j^^r the central density, Mb is the total mass of the bar and 

x 2 y 2 z 2 

m 2 = -^- + 7^ + ^-, a> b> c>0, (6) 
a L b L c L 

with a,b and c are the semi-axes and Mg the mass of the bar component. The 
corresponding potential Vg and the forces are given in D. Pfenniger (1984) ll27ll . 



Applications in the 2D case 

We first apply the SALI index to the 2 degrees of freedom barred potential. In this 
case, we can have a Poincare Surface of Section (PSS) and can check the effectiveness 
of the SALI, by comparing the results. 

In figure 1 , we present two orbits of different kind: (i) (x, y, p x ,Py) = (1.5,0, p x {H) , 0) 
and (ii) (x,y,p x ,p y ) = (— 0.9,0, p x (H),0). For this application H = —0.3. In panels a) 
and b) of figure 1, we show the orbit projections in the (x,y) -plane and in panel c) we 




draw their corresponding PSS, where we can see that the chaotic orbit (i) tends to fill 
with scattered points the available part of the (x,x) plane and that the ordered orbit (ii) 
creates a closed invariant curve. Finally, in panel d) we apply the SALI method for these 
two orbits. For the chaotic one, the SALI tends to zero (10~ 16 ) exponentially after some 
time steps while for the regular orbit, it fluctuates around a positive number. By choosing 
initial conditions on the line x = of the PSS and calculating the values of the SALI, 
we can detect very small regions of stability that can not be visualized easily by the PSS 
method. Repeating this for many values of the energy, we are able to follow the change 
of the fraction of chaotic and ordered orbits in the phase space as the energy of the orbit 
varies. 



Applications in the 3D case 

In figure 2 we present percentages of chaotic, intermediate and regular trajectories, 
where we vary the mass of the bar component (panels A3-B3) and the length of the 
short z-axis (panels A2-B2) of the initial models Al and Bl. The two rows differ in 
the way we give the 27000 initial conditions. For the Al, we give initial conditions in 
the plane (x,p y ,z) with (y,p x ,Pz) = (0,0,0) and for the Bl, in the plane (x,p y ,p z ) with 
(yiZ,p x ) = (0,0,0). By comparing the results, we see that the increase of the bar mass 



causes more chaotic behavior in both cases (panels A3, B3). This confirms the results by 
Athanassoula et al. (1983) 12811 in 2 degrees of freedom. On the other hand, it is obvious 
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FIGURE 2. Percentages of regular (first and second bar) and chaotic (third and forth bar) orbits, for our 
standard model (first row of panels), a model with a thick bar (second row) and a model with a massive 
bar (third row). The two columns show two different ways of choosing the orbital population. 
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FIGURE 3. The percentages of the chaotic and ordered orbits for the two different values of the pattern 
speeds. Panel A corresponds to £2/, = 0.0367032 and panel B to Q,/, = 0.0554349. The latter presents more 
chaotic behavior. 



that when the bar is thicker, i.e. the length of the z-axis larger, the system gets more 
regular. 

Finally, we calculated the same percentages for different pattern speeds Q./,. From 
the orientation of periodic orbits, Contopoulos (1980) Il29ll showed that bars have to 
end before corotation, i.e. that ri > a, where ri the Lagrangian, or corotation, radius. 
Comparing the shape of the observed dust lanes along the leading edges of bars to that 
of the shock loci in hydrodynamic simulations of gas flow in barred galaxy potentials, 
Athanassoula (1992a,b) [ 30],[ 3 1 ] was able to set both a lower and an upper limit to 
corotation radius, namely = (1.2±0.2)a. This restricts the range of possible values 
of the pattern speed, i.e. Q.^ = 0.0367032, that corresponds to the Lagrangian radius 
ri = \Aa and Q.b = 0.0554349, that corresponds to ri = 1.0a. Using the extremes of 
this range, we investigated how the pattern speed of the bar affects the system and found 
that the percentage of the regular orbits is greater in slow bars. 



CONCLUSIONS 

In this paper, we applied the SALI method in the Ferrers barred galaxy models of 2 
and 3 degrees of freedom. We presented and discussed our results comparing the SALI 
index with traditional methods, such as the PSS method for the 2 degrees of freedom and 
showed its effectiveness. We also, calculated percentages of chaotic and regular orbits 
and how they change with the main model parameters, for the 3 degrees of freedom case. 
In more detail, we calculated percentages of chaotic and ordered trajectories, as some 
important parameters, as the mass, the length of the short z-axis and the pattern speed of 
the bar vary in the initial basic model. 
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